Parametric internal waves in a compressible fluid 
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We describe the effect of vibration on a confined volume of fluid which is density stratifled due to its 
compressibility. We show that internal gravity-acoustic waves can be parametrically destabilized by 
the vibration. The resulting instability is similar to the classic Faraday instability of surface waves, 
albeit modified by the compressible nature of the fluid. It may be possible to observe experimentally 
near a gas-liquid critical point. 

PACS numbers: 43.28.+h, 92.60.hh 



It is well known that when an incompressible, viscous 
fluid with a free surface undergoes vertical oscillations, 
standing waves can be parametrically excited which oscil- 
late at half the frequency of the forcing [T] . This Faraday 
instability occurs above a nonzero threshold acceleration 
which is related to the viscosity of the fluid [5]. The 
patterns produced by the interactions of the waves are 
much-studied example of nonlinear pattern formation [3] . 
Driven-dissipative parametric waves are described generi- 
cally by a damped Mathieu equation [2l|4]. In this paper, 
we describe an application of this very general idea to in- 
ternal waves in compressible, density stratified fluids. We 
find that such a fluid can naturally support parametric 
waves and that such waves may be experimentally ob- 
servable near a gas-liquid critical point. This scenario has 
some important differences from previously studies cases 
of parametrically driven internal waves in incompressible 
fluids [SJ |5J |7] because different modes of motion and 
different restoring forces arise in the compressible case. 

The literature on internal waves in density stratified 
fluids is vast, owing to their relevance to atmospheric 
and oceanic processes |8j. Gravity waves, in which the 
restoring forces are due to the weight of the displaced 
fluid, propagate at a flxed angle to the direction of grav- 
ity which is determined by the wave frequency [S]. Ex- 
perimentally, they may be parametrically excited by ver- 
tically vibrating the container [7j, and sidewall reflections 
can give rise to interesting self-reinforcing attractors [6] . 
In these experiments, the density stratification was pro- 
duced compositionally, using dissolved salt. Density gra- 
dients due to temperature gradients can also be excited 
parametrically by shaking the container, although the 
resulting waves may be mixed with other flows due to 
buoyancy driven convection [5] . These previous studies of 
parametrically driven internal waves considered approx- 
imately incompressible fluids, that is, fluids in which the 
flow velocity associated with the wave motion was as- 
sumed to be negligible compared to the speed of sound. 
Here, we allow for the compressibility of the fluid, and 
thus allow for waves restored by non-hydrostatic pressure 



gradients, i.e. sound modes, in the fluid. 

Even nominally incompressible fluids, like water, may 
experience tiny flows under vibration, due to their non- 
zero compressibility. Numerous studies [101 EJ HH HSJ 
[HI [TSl [TB] of shaken fluids have been undertaken both 
on Earth and in space, motivated by the non-ideal mi- 
crogravity environment on orbiting spacecraft. These ef- 
fects, known collectively as g-jitter, are a significant bar- 
rier to using space-based platforms to study ideal gravity- 
free processes. The experiments of Tryggvason et. al. |16j 
showed small coherent motion of neutrally buoyant parti- 
cles in water under vibration in a sealed, isothermal cell. 
These motions have the right scale to be caused by resid- 
ual compressibility effects, and are difficult to explain 
any other way. Such tiny flows are normally negligible 
in experiments, but were here seen superposed on the 
brownian motion and diffusion of microscopic particles 
in a fluid that would be motionless if it were perfectly 
incompressible. Such flows could have an important in- 
fluence on small scale mixing. 

We have theoretically investigated a slab of compress- 
ible fluid subjected to both steady and oscillatory com- 
ponents of imposed acceleration. We show that the os- 
cillatory component can parametrically drive acoustic- 
gravity waves which are supported by the density strat- 
ification due to the constant component of the acceler- 
ation. Our treatment generalizes a number of previous 
calculations [17l [HI [19] , which considered various limit- 
ing cases. 

BASIC EQUATIONS AND SCALING 

In this section we set up the hydrodynamic equations 
that describe the problem. We then consider the appro- 
priate scaling of these equations for the rather unusual 
situation of a fluid which is both compressible and vis- 
cous. The system consists of a laterally infinite layer of 
fluid, confined between two parallel surfaces on which we 
assume rigid boundary conditions. The thickness of the 
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layer is d. The fluid is isentropic and compressible, with 
speed of sound c and molecular viscosity rj. The layer is 
subject to linear accelerations in the direction normal to 
the surfaces. For simplicity, we assume the combination 
of a constant and a sinusoidal acceleration, but the the- 
ory could easily be generalized to more complex linear 
accelerations. 

Conservation of mass requires |30] 



Vu 



IDp 



(1) 



where D/Dt ~ d/dt + {u ■ V) is the material derivative, 
p is the density of the fluid and u = (u, v, w) is the fluid 
velocity. 

Momentum conservation is described by the Navier- 
Stokes equations: 



Du 



P 



Dt 



j77V(V-u)-pz [oo 



aiCOs(a;t)] . 
(2) 

Here, z is a unit vector normal to the surfaces, ao is a 
constant acceleration, while ai is the amplitude of the 
sinusoidal driving acceleration. On Earth, ao = g, and z 
points upward. 

The compressibility of the fluid implies that the pres- 
sure p and the density p are related by an equation of 
state such that 
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where c is the speed of sound and the subscript S refers 
to constant entropy. To recover the incompressible limit, 
the speed of sound c — > oo and Dp/Dt — > 0, such that 
their product remains flnite. The appropriate boundary 
conditions on the velocity are the "no slip" conditions, 
u(z = ±d/2) = 0. 

In order to recast the problem in dimensionless form, 
we must choose scalings for the various quantities in 
Eq. (T]|3l. Under physically reasonable vibration ampli- 
tudes, we expect the flows to be very weak, and therefore 
that the Reynolds number TZ will be small. On the other 
hand, compressibility implies that the Mach number 
while also very small, must be nonzero. Accordingly, we 
choose a scaling such that the pressure gradient term bal- 
ances the viscous term in Eq. ([2]), as would be expected 
for Stokes' flow. We are lead to the following set of di- 
mensionless variables, denoted by primes: 



{x\y\z') 


= {x/d,y/d,z/d), 


t' 


— ty^LUc/d, 


u' 


= u/^/{ujdc), 


P' 


= pVd/{j]^/ujc), 


P' 


= pc^Vd/{r]^/ujc), 


a'o 


= ao/(wc). 


a[ 


= ai/(wc). 



(4) 



In terms of these variables, omitting the primes, Eq. ([ij 
is unchanged, and the other basic equations are 
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(5) 
(6) 



Dp _ Dp 
iDt ^ Ijt' 

where the Mach number M — tod/c. The Reynolds num- 
ber under the same scaling is simply TZ ~ lodjc = M.. 
Thus, it is possible for both M. and TZ to be small, while 
their ratio is 0(1). This is a seldom used limit in fluid 
mechanics: that of Stokes' flow at nonzero Mach num- 
ber. However, Eqs. ([l]) and ([s]), ^ are still completely 
general, and no assumptions have been made about the 
size of M. . 

It is convenient to absorb a factor of M. into the di- 
mensionless acceleration parameters and deflne 

Ao = Alao = dflo/c^, and ^ Ma!^^ da^j (? . (7) 

We will show that these elementary equations contain 
sound modes and internal gravity wave modes which can 
be parametrically excited by the driving acceleration. 

The base state corresponds to the case A\ = 0, with 
Aq arbitrary. In this case, the fluid is static and there is 
a hydrostatic pressure distribution given by 

J = -A.p. (8) 

The corresponding density distribution is given by 

poiz) = Po(0)e-^°^ (9) 

where po(0) is the density at z = 0. The base state 
is thus density stratifled with scale height Aq^ , which 
could be much larger than the thickness of the layer [35 . 
The typical order of magnitude of change in density of 
normal liquid like water in a 10 cm box in the ground 
based experiments is rather small and changes by only 
O(10-4)%. 

Density stratification due to compressibility is funda- 
mentally different than that due to salinity, the situa- 
tion most often used in laboratory studies of internal 
waves Nevertheless, such stratification is expected to 
become significant near the critical point of gases, where 
the speed of sound goes to zero and the compressibility 
diverges [3T1[52]. In practice, the details of the equation 
of state of the gas are required to calculate the stratifica- 
tion. Close to the critical point, the density profile will 
become quite nonlinear |31j . 



LINEARISED EQUATIONS AND NUMERICAL 
SIMULATION 



We now consider the fluid under vibration and derive 
the equations for the growth or decay of small deviations 
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about the "base state" . We consider only a linear theory 
valid in the limit of small deviations. This leads to a 
well-known equation for parametric oscillations. 

To test the stability of the base state we consider 
the fluctuations (5u, 5p and 5p oi the velocity, pressure 
and density fields and linearise the equations of motion, 
Eqs. Q, ([5]) and Q. Linearisation of Eq. ([!]) about the 
base state with p = po and u=0 gives 



poV • Su - 



dSp 
~dt 



Sw 



dpo 
dz 



0. 



(10) 



Considering the liquid to be irrotationaL, i.e., applying 
V(V • 5u) — V^Su, linearisation of Eq. ([5| results in 



Mpo 



dSu 

~dt' 



= -V5p+ -W'^Su 

- Sp[Ao + Aicos{y/Mt)]z 

— pqAi cos{VMt)z. 



(11) 



Finally, the equation of state reduces to 
dSp 85 p 



dt 



dt 



(12) 



It is straightforward to eliminate Sp and dp between 
Eqs. (10 1, (111 and (12 1. After some simplification and 



assuming Ai and A4 are small so that the inhomogeneous 
forcing term can be neglected we obtain 



M 



d^Su 1 

Air 



, dSu 

Ik 



+ AoV(5 • Su) = V^(5u 



+z- 



Pa 



PoV ■ Su + {Su ■ V)po cos{VMt), (13) 



If Aq = 0, so that the constant component of the accel- 



eration is zero, Eq. ( 13 1 describes damped acoustic waves. 



It can be shown from the dispersion relation for inviscid 
fluids that Su and the wave vector k are in the same di- 
rection and thus that such waves are purely longitudinal. 
In the opposite limit, when Ai — and the oscillatory 



component of the acceleration is zero, Eq. ( 13 ) reduces 



to a dynamical equation that supports mixed acoustic- 
gravity wave motion which is damped by viscosity. Such 
mixed waves are no longer completely longitudinal |30j . 



Eq.(13l was studied in two limited cases by |T7],[T5] 
and [T5]. In [T7j and [TS] the effects of compressibility 

on the early stages of Rayleigh- Taylor instability in non- 
viscous liquids was discussed. In [I9j it was shown that 
in absence of gravity for many liquids with low viscos- 
ity (e.g. air, water) Eq. (13 1 can be approximated by a 
damped wave equation. Moreover in the following anal- 
ysis it is shown that this damping term for low viscous 
liquids has negligible effect on the stability of the strat- 
ified state in presence of jitter. However for brevity we 
kept this term in our analysis. 

Eq.(13l is a parabolic, third-order partial differential 
equation. Homogeneous form of Eq.(13) without jitter 



for harmonic waves are solved by Stokes (5U] Stefan [5T] 
and Rayleigh. Stokes transient problem has two types 
of solution : (a) solutions in the form of series [22l [23] 
or integrals [21] and (b) closed-form approximations [2S] . 
Cobbold et al. [23] has recently stated that a difficulty 
arises with the approximate, closed-form type of solu- 
tion : " approximate solutions to such problems do not 
satisfy causality in the strict sense ." Many examples of 
breakdown of causality for transient solutions of Stokes 
equation may be found in White pT. Jordan et.al [TS] 
and Blackstock|25j have raised doubts over the validity of 
Stokes equation itself [2^], with Jordan et. al [T^ stating 
that solutions for this problem "...do not satisfy causal- 
ity" . They along with others [23j [25] have also observed 
that transient solutions of Stokes equation "are felt in- 
stantly" throughout the entire ffuid domain. It has been 
established that an instantaneous response throughout 
the fluid is unphysical, since an infinite wave speed is not 
permitted. Weyman [28], pointed out that this not only 
Stokes' equation but also other parabolic equations such 
as, the diffusion equation, is commonly said to lead to an 
infinite speed of propagation through the medium sup- 
porting the motion, even though such behavior is entirely 
unphysical. Thus, it seems that the parabolic differen- 
tial equations, even though derived correctly may lead to 
unphysical results. To avoid this paradox we followed 
Jordan et. al [TS] and approximated the damping term 
to convert the parabolic equation into a well-posed hy- 
perbolic equation. For small Aq and Ai we assume the 
approximation |19j is still valid and we have approxi- 



mated Eq.( 13 ) as 



M 



d^Su dSu 



df^ 3po dt 
.Al 



AoV{z ■ Su) 



Po 



poV • Su + {Su ■ V)/9o cos{V Mt), (14) 



It follows that for liquids bounded by two parallel plates 
a solution of the form Su ~ u{t)exp {az -f i{kix + k2y)) 
can be sought. Substituting this form and restricting our- 
selves to the long wavelength modes (fci, ^2 0) we get 
an equation for the velocity component in the z direction 
i.e., Sii as 



M 



d'^w{t) ^ dw{t) 

dt^ 3po dt 
-Ai{a ~ Ao)cos{VMt) w{t) 



a{AQ — a) 
0, 



(15) 



where w{t) is the time varying part of Sw. This equa- 
tion is a modified form of the Mathieu equation which 
describes parametric oscillations. It is evident that a 
parametric response to the driving frequency is possi- 
ble for a range of parameters M, Aq, Ai and po, for a 



nonzero value of a. We have solved Eq. (15 1 numerically 



using a four point Runge-Kutta method and the results 
are shown in Fig. [l] We find that in presence of constant 
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acceleration, and in absence of shaking {Aq 0, Ai = 0) 
the stratified fluid can undergo internal oscillations at the 
characteristic frequency N, where N = ^ a^A^ — at). 

The amplitude of the response is modulated by the 
presence of shaking {A\ ^ 0). The system becomes para- 
metrically unstable if the amplitude of shaking exceeds 
that specified by the stability boundary derived analyti- 
cally in the next section. 



PERTURBATION ANALYSIS 



Eq. (15) with A\ = 0, i.e., in the absence of shaking, 



will always have oscillatory solutions for Aq > a. But any 
viscosity, however small, will obviously cause the oscilla- 
tions to damp out because the coefficient of the second 
term in Eq. ( 15 1 is always positive. In what follows 



we 



will show that such an oscillation of the density strati- 
fied fluid can grow under the action of periodic forcing 
or shaking. 



Let us first rewrite Eq. ( 15 1 as 



+ if^l^ + (/? + e cos(V7Wi)) = 0, (16) 



where K = M/ipo, (3 — a{Ao — a)/M and e = 
Ai{AQ — a)/M. For small e, it is possible to use 
very standard methods to do a perturbation analysis of 
Eq. ( [T6| . This analysis finds the onset of linear instabil- 
ity and the regions where solutions of Eq. (16 1 become 



unbounded in the linearized theory. The only restriction 
imposed on the parameters is that |e| ^ 1. For small |e|, 
we expand (3 and Sw in powers of e as 



w 



wq + ewi 



and 



(17) 



(18) 



lowing equations at 0(1), 0(e) and 0(?): 



Putting the above expansions in Eq.(|16|) we get the fol- 

(19) 



dt^ 



at 



d?wi dwi 



d^W2 

df^ 



^dW2 ^^^^ ^ -P2W0 - PiWi 

at 



■ cos{VMt)wi 



(20) 



(21) 



The solution of Eq. (19 1 is of the form wq = b^e 



At 



with A = —K 12^^ K'^ — 4/?o/2, where 60 is an arbitrary 




FIG. 1; Response of fluid particles with time, a) periodic 
response of fluid particles when density stratified base state 
due to compressibility is perturbed in absence of jitter. A\ = 
0, Ao = 1.0, = 0.3. b) amplitude is modulated by the para- 
metric forcing term. A-y = 0.03, = 1.0, a = 0.1, = 0.3. 
c) response for A-y = 0.05, Aq = 1.0, a = 0.1, = 0.3, dis- 
placement of fluid particles remains bounded, d) solution be- 
comes unbounded when the value of A\ is above the stabiUty 
boundary specifled by Eq. ( 26 1 . Linear theory breaks down in 
this case. 



constant. For low viscosity liquids, K is small and when 
< 4(3q there will be an oscillatory component in the 
leading order solution. 

When Wq is put into Eq. (20 I, the 0(e) equation, the 



first term on the right hand side contains the secular 
term and has to be zero for the 0(e) solution to exist, 
according to Fredholm's condition Thus, /3i — 0, 

and the solution at 0(e) is 



wi = b° cos{VMt)e^^ + hi sm{VMt)e 



.At 



(22) 



where = bo /[I + {K + 2Xy 



and bl 



iK + 2X)b°. 



Now, putting Wq and Wi into Eq. ( |20| and removing 
the secular term we find P2 = -(^1/2) and /3 calculated 
up to this order is 



f3 = f3o 



bo 



2[1 + {K + 2XY 



(23) 



The O(e^) solution shows that the amplitude of the peri- 
odic response will be modulated by the parametric forc- 
ing term. 

Next, we consider the transformation 



wit) 



w'it) 



(24) 



in order to cast the damped Mathieu equation, Eq. ([16|, 
into the form of a normal undamped Mathieu equation. 
The above transformation means that only those un- 
bounded solutions of the Mathieu equation which grow 
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faster than e~^*/^ will actually destabilize the oscillatory 
stratified state. An analysis along standard lines using 
Floquet theory [34] reveals the first unstable region of 
the stratified state to occur when 



(/3- 



M 



1 



y - -{e' -MK^) < 0. 



(25) 



From the above inequality, the minimum shaking ampli- 
tude to induce unbounded motion turns out to be 



Al > 



A[a{Ao -a)- (MV4)]- 
{a-Aof 



(26) 




FIG. 2: The stability regions in the forcing amplitude Ai 
vs. Mach number M plane, for a fluid with negligible viscos- 
ity. The region above the curves is parametrically unstable, 
and the internal waves will grow to become nonlinear. The 
sharp minimum means that forcing with infinitesimal ampli- 
tude may destabilize the system at those parameter regimes. 



also interesting that waves can be driven parametrically 
with modest values of Ai, even at frequencies well below 
N/2. 

It is important to note that standard results for the 
damped Mathieu equation indicate that higher order 
resonances will be encountered for higher positive values 
of the quantity cr = /? — (i^^/4), where 



a{Ao - a) M'^ 



M 



(28) 



The second term on the r.h.s. of Eq. (28 1 is negligible 



for low viscosity fluids. General stability diagrams which 
show the higher order resonances in the a ~ K plane are 
beyond the reach of perturbation theory and must be ob- 
tained numerically [S^. These higher order regions have 
considerable spread in the parameter space, and thus for 
a higher forcing amplitude it is more probable that the 
system will be in a region of unbounded solutions. 

(28), the first term is positive when Aq > a 



In Eq. 



and can slowly be increased if M. is gradually decreased. 
This suggests a possible experimental protocol to fix the 
driving frequency and acceleration and slowly increase 
the temperature of the gas away from near the gas-liquid 
critical point, where the speed of sound becomes small. 
This has the effect of decreasing M. One could then 
in principle successively observe the various higher or- 
der resonances which are characteristic of the Mathieu 
equation. A detailed calculation using real gas equations 
of state, combined with numerical solutions of Eq. (15 1 



would be necessary to make specific predictions for ex- 
periments, however. 



CONCLUSION 



We can rewrite this result in dimensional form as 



4c4[a(ao-ac')-(c^V4)]^ 
(ac2 -aof 



(27) 

where a = (x/d and k = k/d. For many liquids like 
water the first term on the r.h.s. of Eq. (27 1 is negligible 
since kinematic viscosity v = rj/pQ ^ 10^'^ m/s is quite 
small. The second term then dominates and determines 
the minimum forcing amplitude. The corresponding forc- 
ing frequency is given by ujq = 4a(ao — ac^). Typically, 
there will exist a narrow band of frequencies for which 
the system is unstable just above this minimum forcing 
amplitude. 

In Fig.|2] we plot the stability threshold for Ai against 
M for different values of Aq, assuming a fluid of negligi- 
ble kinematic viscosity. We find that driving frequencies 
close to a subharmonic resonance with the characteristic 
frequency N most easily drive the internal waves. It is 



We have investigated the effect of vibration on a com- 
pressible fluid which may also be subject to a constant 
acceleration. The response of the fluid can show both 
harmonic and subharmonic components, as well as pos- 
sible higher order resonances. The harmonic response 
probably explains the coherent oscillations observed by 
Tryggvasson et al [16] in a sealed water cell. The os- 
cillatory acceleration can interact with the stratification 
produced by the constant acceleration to produce para- 
metrically excited internal waves. For fluids of small vis- 
cosity, the threshold for this instability can be very small 
if the frequency is selected to be close to a subharmonic 
of the Brunt- Vaisala frequency of the stratification. All 
these effects may be observable in fluids near their critical 
points, where compressible effect become very large [55] . 
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